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Semiconductor devices have been scaled to the point that transport can be dominated by only 
a single dopant atom. As a result, in a Si Fin Field Effect Transistor Kondo physics can govern 
transport when one electron is bound to the single dopant. Orbital (valley) degrees of freedom, 
apart from the standard spin, strongly modify the Kondo effect in such systems. Owing to the 
small size and the s-like orbital symmetry of the ground state of the dopant, these orbital degrees 
of freedom do not couple to external magnetic fields which allows to tune the symmetry of the 
Kondo effect. Here we study this tunable Kondo effect and demonstrate experimentally a symmetry 
crossover from a SU(4) ground state to a pure orbital SU(2) ground state as a function of magnetic 
field. Our claim is supported by theoretical calculations that unambiguously show that the SU(2) 
symmetric case corresponds to a pure valley Kondo effect of fully polarized electrons. 

PACS numbers: 71.27.+a, 71.30.+h, 73.23.Hk,72.15.Qm 

The resistance of metals with magnetic impurities anomalously increases as one decreases the temperature. This 
Kondo effect pQ can be explained as the screening of the localized spin of the magnetic impurity by the spins of the 
de-localized electrons in the metal. As a consequence of this screening, the localized spin and the itinerant ones form a 
many-body singlet with binding energy Tr-, which defines the low temperature scale at which Kondo physics appears. 
A few years ago, it was shown that quantum dots (QDs) [2] behave as Kondo impurities. The transport properties 
of QDs in the Kondo regime are quite remarkable: starting from an insulating QD in the Coulomb blockade regime 
at high temperatures, the linear conductance reaches the maximum unitary value of a perfect quantum conductor, 
namely G = dl / dVb\v b ^o — 2e 2 //i as the temperature is reduced well below Tk OP- At finite bias voltages V&, Kondo 
physics manifests as a zero-bias anomaly in the dl/dV^ curves whose width is roughly given by Tk- The Kondo effect 
in QDs originates from quantum fluctuations of the charge residing in the QD: electrons can transit through virtual 
states on a time-scale which is shorter than allowed by the Heisenberg uncertainty principle [1 j. This mechanism 
generates effective spin flips which in turn lead to Kondo physics. Importantly, the role of the electron spin can be 
replaced by any other quantum degree of freedom such as e.g. orbital momentum [3H2], giving rise to exotic Kondo 
effects. Furthermore, the simultaneous presence of both a spin- and an orbital- degeneracy leads to an SU(4)- Kondo 
effect, where SU(4) refers to the symmetry of the corresponding Kondo ground state [3H9]. 

In the past, SU(4) Kondo symmetry has been predicted to arise in parallel double quantum dot systems [3], but 
so far it has only been clearly observed in carbon nanotubes [5] and in single dopant devices in Si [9]. Si is a good 
candidate for observing SU(4) Kondo physics due to its six- fold valley (orbital) degeneracy of the conduction band 
and orbital effects are the most probable cause of the unexplained (at the time) behaviours observed in early studies of 
Kondo in Si QDs [TO]. Here we show a fully tunable Kondo effect in a Si Fin Field Effect Transistors (FinFETs) [TIHT4] 
[see inset of Fig. 1(a)]. We use a gate voltage to bring the orbital states into degeneracy, while independently tuning 
the spin splitting by means of an external magnetic field. Interestingly, in our system, orbital Kondo physics survives 
at very high magnetic fields (even for Az > ^bTk^ with Az = g*/isB being the spin splitting) allowing us to tune 
our device from an SU(4) to an SU(2) symmetry. In this latter case, our results can be understood as a pure orbital 
Kondo effect of spin-polarized electrons [4] [15] . Our claim is fully supported by theoretical calculations that take into 
account both thermal and quantum fluctuations in a non-pertubative way. 

Single dopants can be individually addressable in FinFETs [TTHT4] . Our devices consist of Si nanowire connected to 
large contacts etched in a 60 nm layer of p-type Silicon On Insulator (SOI). The wire is covered with nitride oxide (1.4 
nm equivalent Si02 thickness) and a narrow poly-crystalline silicon wire is deposited perpendicularly on top to form 
a gate on three faces. Doping by ion implantation with As over the entire surface forms n-type degenerate source, 
drain and gate electrodes while the channel protected by the gate remains p-type. The conventional operation of this 
n-p-n field effect transistor is to apply a positive gate voltage to create an inversion in the channel and allow a current 
to flow. Unintentional As donors may be present below the Si/insulator interface showing up in the sub-threshold 
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FIG. 1: a) Stability diagram of the differential conductance (for B — T and T— 290 mK) of a single atom transistor. The 
transport resonances inside the diamond, which are Kondo related, should not be allowed due to Coulomb blockade (CB). In 
the inset, a Scanning Electron Micrograph (SEM) of our device with the schematic position of the dopant in the channel (in 
red), is shown, b) Enlargement of the Kondo region of Fig. Q](a). c) Evolution under magnetic field of the Kondo resonance 
at 290 mK and for V g = 464 mV (far from the degeneracy point). For this gate voltage value, the expected 2g/i J eB splitting 
of the central Kondo resonance is observed and as expected for non-degenerate SU(4) Kondo, the side resonance do not split 
in magnetic field 9 . d) Same as Fig. [I] (b) but for B— 10 T. In this figure it is possible to observe an unusual (compared to 
other Kondo regions) survival of the zero bias peak in the V g ~ 450 mV region. 



transport characteristics [TT| H2]. 

From a large ensemble of FinFET devices, we select those that show the fingerprint of single donor transport, which 
essentially consists of a pair of characteristic transport resonances associated with the one-electron (D°) and two- 
electron (D~) charge states of the single donor [TTHT3] . The electron ground state originates from the hybridisation 
of the donor hydrogen-like state with a quantum well state formed at the Si/insulator interface by the high electric 
field in the channel [TTJ [12j [161 Hi] • Valley degeneracy in Si is strongly broken at the donor (valley splitting ~ 20 meV 
[18] ) but a nearly two- fold degeneracy remains at the interface. As a result, the two lower orbitals show a splitting 
A typically on the order of a few meV [12]. This valley splitting may be modified externally by applying a gate 
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FIG. 2: Current measured in the CB region of Fig. [T] (a) for zero bias voltage, T — 288 mK and increasing values of B. The 
peak observed for B = 10 T and around V g ~ 450 mV (red circle) is related to pure orbital fluctuations while the red line 
indicates the crossover from a Kondo resonance with SU(4) symmetry and degenerate levels to a pure SU(2) orbital Kondo. 
The Kondo temperatures for B = T and V g = 450 mV, 455 mV and 464 mV have been estimated to be ~ 8 K, 6 K and 4 
K respectively. An increasing offset has been added to the curves for clarity. 



voltage [12] which affects interface quantum well depth and the transparency of the barriers leading to a modification 
of the energy levels of the hybridized electron ground state wave- function [19 , 20 . On the other hand, the s-like orbital 
degree of freedom is still dominant in this system [12], and, consequently, the valley splitting is virtually independent 
to magnetic fields applied parallel to the channel. 

Transport spectroscopy measurements as a function of gate and bias voltage have been performed using a lock-in 
technique, typically with a 50 /iV amplitude ac voltage, and a 3 He cryostat with a base temperature of 270 mK. A 
magnetic field B parallel to the direction of the channel of the FinFETs has also been applied. All the data presented 
here correspond to a device in which a single dopant atom can be individually addressed due to strong capacitive 
coupling with the source and drain contacts [11] (see Fig. [5] in the supplementary section). 

In Fig. [l](a) the differential conductance G— dl/dVb of our device as a function of the gate V g and bias V& voltage 
is shown at T— 290 mK and B = T. A Coulomb diamond is defined between the D° peak at V g — 426 mV and the 
D~ peak at V g ~ 480 mV [TTJ [12]. As it has been already shown in a previous publication [9 , for our system the 
Kondo resonances around Vb ~ have the typical SU(4) Kondo signature. In particular, we observe a split Kondo 
resonance associated to fluctuations in two quasi-degenerate orbital states and a central one associated with spin 
fluctuations. This is illustrated in the blow up in Fig. [l](b) which reveals a non-zero bias central peak and side peaks, 
all of these scaling in temperature according to the same empirical law [21] , as described in Ref. [9] . The spin-related 
peak occurs at non-zero bias due to asymmetries in the contacts [20 , 22 (see also theory section in the supplementary 
section) and to the possible presence of Fano physics [23j [24] , the signature of which has been already observed in our 
systems [25], while the side peak arises due to the valley splitting (A ^ mV) of the donor bound electron ground 
state. The evolution of the Kondo peaks with an increasing magnetic field for V g ~ 464 mV is shown in Fig. [I] (c). 
The strongest peak splits 2g* [IbB, as expected [2], while the satellite valley peak remains unaffected by the magnetic 
field [9]. These features are common for most values of the gate voltage [9 . The behaviour with magnetic field is 
dramatically different for V g ~ 450 mV (see Fig. [I] (d), Fig. [5] and Fig. [3] (a)). This gate voltage is special as it 
corresponds to a valley degeneracy point (A = mV) [9]. At gate voltages V g = 450 mV, 455 mV and 464 mV, the 
splitting of the orbit als progressively increases, which results in a decrease of the corresponding Kondo temperature 
[3] [26], i.e:. Tk = 8 K, 6 K and 4 K respectively. Therefore, for V g ~ 450 mV, we have a four- fold degeneracy that 
reflects in a SU(4) Kondo effect. In the absence of a magnetic field, the dominating resonance is a structured peak 
at non-zero bias: Vb = -1 mV. However, as a magnetic field is applied, this non-zero bias resonance suffers a 2g/i#B 
splitting and gets suppressed while a simple central (Vb ^ mV) peak gets enhanced and dominates at large fields 
B > 4 T. These behaviours are depicted in Fig. [3] (a), with the differential conductance as a function of Vb and the 
magnetic field, and in Fig. [3] (b) with the evolution of the heights of the two peaks, at Vb — — 1 mV and Vb — mV, 
as a function of the magnetic field. From these data we can observe that, for B = 10 T, a central resonance (with 
Tk ~ 6 K, see Fig. 5 of supplementary section) dominates the stability diagram (as also shown in Fig. [l](d)). The 
development of this central (Vb — 0) peak with increasing magnetic field is also illustrated in Fig. [2j where the current 
versus V g is plotted for different values of B. As the black arrow in Fig. [3] (a) indicates, this peak does not shift or 
split in Vb when increasingly high values of B are applied, implying that it is related to the orbital (valley) degree of 
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FIG. 3: a) Evolution of the Kondo peak at 288 mK and for V g = 450 mV. For this gate voltage value no splitting is observed 
for B > 4 T as the central peak survive (black arrow) also for very high magnetic fields, b) Traces that show the evolution 
under increasing magnetic field of the peaks height at Vb ~ -1 mV (rhombuses) and of the peak at V& = mV (circles). 



freedom which must be preserved during tunnelling [6j |9] . 

This observation is of importance as the conservation of valley index upon tunnelling is most probably sample 
dependent and still the subject of theoretical debate HUES- Moreover, the large value of the magnetic field indicates 
the spin must be completely polarized [U [15] , such that the only Kondo active degree of freedom is the orbital one 
(implying an SU(2) symmetry as schematically shown in Fig. [4](c)). Last, it is possible to interpret the white circles 
m Fig. H (b) as G versus T/Tk (as B increases, Tk reduces so T/Tk increases). The observed G versus T/Tk 
dependence indeed follows the non-monotonic behaviour expected for a two-stage Kondo system, in agreement with 
the theoretical description of the SU(4) to SU(2) symmetry crossover [6| Yf\ 126]. 

Having discussed the main experimental features, we now focus on a theoretical analysis that supports our previous 
interpretations. Our theoretical calculations are performed by considering a QD with two (nearly) degenerate localized 
orbit als (valleys) coupled to reservoirs. The effect of the external magnetic field is to lift the spin degeneracy of the 
single-particle energy levels, with the resulting splitting denoted Azeeman (see Fig. [4] (a) and [4] (b)). The final 
results are obtained by using the so-called non-crossing approximation (NCA), a fully non-perturbative theory which 
includes both thermal and quantum fluctuations, complemented with the Keldysh-Green's function method to take 
into account non-equilibrium effects, due to the finite bias voltage Vb. More details on the theory can be found in the 
supplementary section. In Fig. [4] (a) and|4](b), the calculated differential conductance as a function of Vb and the 
magnetic field is plotted. At zero magnetic field, we obtain a SU(4) Kondo resonance which is maximum at non-zero 
bias voltage. This effect is attributed to an asymmetry in the way the dopant is coupled to the reservoir [20 j. We treat 
the case of orbital degeneracy, namely A = mV. As a magnetic field is applied, this non-zero bias resonance splits 
in three: the two outer peaks correspond to inelastic cotunnelling-like processes in which spin flips occur mediated 
by Kondo fluctuations in the orbital degree of freedom. The central peak corresponds to orbital Kondo processes in 
which the spin is conserved. As the magnetic field increases, the central peak moves towards zero bias and the outer 
peaks have no longer the same height (the one corresponding to positive voltages being higher). This, again, can be 
easily understood as resulting from coupling asymmetry: for positive voltages, the distance between the spin-polarized 
orbital dopant levels and the left chemical potential decreases resulting in a stronger Kondo effect as compared with 
the one for negative voltages where the right reservoir (the one with lower chemical potential) is less coupled to the 
dopant levels. Importantly, only quantum fluctuations between orbital levels are active so the underlying symmetry 
of the Kondo effect is SU(2). We have assumed that the valley index is preserved during tunnelling. This is a 
necessary condition for the observation of the SU(2) orbital Kondo effect as valley index mixing would destroy it 
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FIG. 4: a) Theoretical investigation of the degenerate SU(4) to pure orbital SU(2) crossover as schematically described in Fig. 
[4](c), performed using the non-crossing approximation (see supplementary section and Ref. b) G versus bias traces 

for increasing values of B along the crossover as in Fig. [4] (a). The black, red, blue, violet, green and yellow curves represent 
respectively the Azeeman = Tk, 10 Tk, 20 Tk, 30 Tk, 40 Tk and 50 Tk cases, c) Schematic of the phenomenom discussed 
in this paper. Starting from degenerate valleys (represented in black and in red in this figure), an increasing magnetic field B 
saturates the spin degree of freedom due to the Zeeman splitting, producing a crossover from an SU(4) Kondo symmetry to a 
SU(2) pure orbital Kondo symmetry. 



[8 . It is important to mention that, in order to obtain pure orbital SU(2) physics in the theoretical calculations, we 
need to include higher magnetic fields as compared to the experiment. This can be easily understood because spin 
decoherence microscopic mechanisms (apart from cotunnelling-like finite lifetimes) are not included in the calculations 
but are surely present in the experiments (see Ref. [9 and supplementary section). To be more specific, finite-bias 
decoherence is also the cause of the experimental suppression of the SU(4)b=ot/ SU(2)b=iot peaks ratio observed 
in Fig. [3](b). This discrepancy with theory, see Fig. [4] (a) and Fig. [4] (b) or Ref. [6 j, has nothing to do with Kondo 
physics but with the fact that the experimental peak at B = 10 T is less affected by decoherence processes. Overall, 
there is a very good agreement between the theoretical results and the experiments, as in both case the survival, even 
at very high magnetic fields, of a simple central (V5 ~ mV) pure orbital Kondo related peak, is observed. This 
supports an interpretation of the experimental data in terms of an SU(4) to SU(2) crossover driven by magnetic field, 
where we tune the system between two distinct Kondo states. This behaviour is also in agreement with experiments 
in carbon nanotube quantum dots in the presence of a perpendicular magnetic field [27] . Similarly to our case, the 
magnetic field in these experiments only couples to the spins and renders the orbital fluctuations unaffected. In 
conclusion, we have demonstrated a controlled crossover between SU(4) and SU(2) Kondo states driven by magnetic 
field in a nanoscale Si transistor. The latter SU(2) Kondo effect originates only from quantum fluctuations in the 
orbital (valley) degree of freedom. We also perform theoretical calculations, including both spin and orbital Kondo 
physics as well as non-equilibrium effects, which confirm our experimental findings. Overall, the importance of the 
orbital degree of freedom in novel nano-scaled Si systems is confirmed, opening the way to the possible use of different 
symmetries that the donor orbitals can provide for innovative implementations of Si quantum-electronics [28] such as 
for example, valley-based quantum computation [29] . 
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ducted within the ARC-CQC2T (project number CE110001027). R.A. and M.J.C. acknowledge funding from MICINN 
(Spain) through Grant No. FIS2009-08744. M.J.C. also acknowledges the Ramon y Cajal program, MICINN (Spain). 
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FIG. 5: Magnetic transport data of the D° and D peaks for the device discussed in the main manuscript. 



I. SUPPLEMENTARY SECTION: MAGNETIC FIELD PROBING OF AN SU(4) KONDO RESONANCE 

IN A SINGLE ATOM TRANSISTOR. 

A. Magnetic Transport Data of the D° and D~ peaks 

In Fig. [5] it is shown that, as expected for a single atom transistor [11], an odd/even effect in the magneto-transport 
spectroscopy data is observed. This consists of the Zeeman down shift of the first electron resonance (D°) and the 
Zeeman up shift of the second electron resonance (D~) [11 . The doubling of the D~ peak is a feature sometimes 
observed in a single dopant system and can be associated with the presence of a nearby defect quickly switching 
its charge state (see Ref. [9] and references therein). These results, together with the charging energy of 30-35 mV 
between the D° and the D~ state in this sample (Fig. [I] (a) of the main text) and the identification of the excited 
states as in Ref. [12 , confirm the single dopant nature of the Coulomb blockade transport observed in the device. 



B. Kondo Temperature for the Central Resonance at T and at 10 T 

In Fig. [6j the evolution with temperature of the central resonance for B = T and V g = 450 mV (centred at 
= mV) and the evolution in temperature for both the B = T and B — 10 T peaks are shown. A value of Tk 
= 6 K is obtained for the B = 10 T case, which is higher if compared to SU(2) spin Kondo temperatures observed 
for quantum dots (typically < IK) [9 j. Furthermore, the point for B = T and T = 4.2 K is considerably lower 
if compared to the values in the 0.3 K-1.6 K range. The peak saturation observed for < IK in the B = T case, 
strongly supports the conclusion that, at T ~ 300 mK, our system is indeed in the zero temperature limit. 
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FIG. 6: a) Evolution with temperature of the central peak at V g ~ 450 mV and for B — 10 T between 0.3 K and 1.5 K. This 
situation is similar to the one theoretically described in Ref. 15 . b) Evolution in temperature of the peaks maximum at B = 
T (in red and with T K « 8 K) and at B = 10 T (in black and with T K « 6 K). To obtain the T K , we have fitted the data 
using the empirical law proposed by Goldhaber-Gordon et al [21]. Instrumental error bars compatibles with a 98 % confident 
band have been used. 



C. Theory Model 

A single dopant in a FinFET can be described as a quantum dot (QD): The single electron state corresponds to the 
neutral donor D° while the two-electron state corresponds to a negatively charged donor D~ . Consider then a QD 
with two (nearly) degenerate localised orbit als (valleys) coupled to reservoirs. Hereafter we will denote this orbital 
quantum number by \i = 1, 2. The dot is then described by the Hamiltonian 

where e^ G is the single-particle energy level of the localized state with orbital \i and spin <j, eft (d^) the fermion 
creation (annihilation) operator of the state, n Mcr = S d^ the occupation, (fi = 1,2) the intra-orbital Coulomb 
interaction, and U12 the inter-orbital Coulomb interaction. The effect of the external magnetic field is to lift the spin 
degeneracy of the single-particle energy levels. The resulting splitting is denoted Az so that the single-particle energy 
levels e^ a have the form e^/i = eo ± (A^/2). 

On the other hand, as we argue in the main text, the valley (orbital) states are s-like [12 and are therefore not 
affected by the magnetic field. The precise values of the Coulomb interactions depend on the details of the 

system, but should be of the order of the charging energy E D - — 2E D o ~ 36 meV [30]. We focus on the regime where 
the system of the localized levels is occupied by a single electron (^^(n^) w 1, quarter filling), and the charging 
energy is much bigger than other energy scales (i.e. U — >• 00). In this regime the Hamiltonian in Eq. suffices to 
describe all relevant physics. 
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FIG. 7: Density of states for Vbi as = with symmetric (left) and asymmetric (right) couplings to the left Yl and right Yr 
reservoirs. 



Kondo physics arises as a result of the interplay between the strong correlation in the dot and the coupling of the 
localized electrons with the itinerant electrons in conduction bands. In our case, these are described as two leads 
(a = L and R) which are treated as non-interacting gases of fermions 



k M=l,2 cr 

When the leads have the same symmetry as the dot, the orbital quantum number /i in the leads is identical to the 
orbital quantum number in the dot and will be preserved over the tunnelling of electrons from the dot to the leads and 
vice versa. This is what we expect to have in our system with heavily doped Si as source and drain and a quantum 
dot that results from the hybridisation of the donor hydrogen-like state with the quantum well state formed in the Si 
channel [12]. This situation is described by the tunnelling Hamiltonian 



H T = ( V *k^al kfI(J d^ + fc.c.) . (3) 

akficr 

The total Hamiltonian is then given by H = Hl+Hr+Ht+Hjj . For simplicity, we ignore the k- and cr-dependence 
of the tunneling amplitudes. Therefore, we consider a simplified model with Va,k^a — Va/V% which defines the widths 
Tq; = 7rpo\V a \ 2 , where po is the density of states in the reservoirs. 



D. Non-crossing approximation method 



Now we write the physical fermionic operator as a combination of a pseudo-fermion and a boson operator as follows: 
d^a = ^ ' /fur where is the pseudo-fermion which annihilates one "occupied state" in the pth localised orbital with 
spin cr, and is a boson operator which creates an "empty state" . Quite generally the intra(inter) Coulomb interaction 
is very large and we can safely take the limit of U oo. This fact enforces the constraint 
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FIG. 8: Density of states for zero Zeeman spliting and different values of Vbi&s with symmetric (left) and asymmetric (right) 
couplings to the left Yl and right Yr reservoirs. A finite bias splits the Kondo peak in 2. 




FIG. 9: Density of states for finite Zeeman splitting and different values of T4ias with symmetric (left) and asymmetric (right) 
couplings to the left Yl and right Yr reservoirs. 
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FIG. 10: Differential conductance for symmetric (left) and asymmetric (right) contacts for different magnetic fields. The 
conductance is given in units of the quantum of conductance e 2 /h (which corresponds to 38. 74/^5, as presented in the main 
text). 



Xl/icr fjiaf^cr + tfb = 1, that prevents the accommodation of two electrons at the same time in either the same 
orbital or different orbitals. This constraint is treated with a Lagrange multiplier. 



#sb= E^ c l,^.-+E £ o.^/^ + -# E (iA'+H^P^^^" 1 • (4) 

Notice that we have rescaled the tunneling amplitudes — » V ^VN according to the spirit of a 1/iV-expansion (TV 
is the total degeneracy of the localized orbital). 

Our next task is to solve this Hamiltonian, which is rather complicated due to the presence of the three operators 
in the tunnelling part and the constrain. In order to do this we employ the so-called Non-Crossing approximation 
(NCA) [3TH33] . Without entering into much detail of the theory, we just mention that the boson fields in Eq. Q 
are treated as fluctuating operators such that both thermal and charge fluctuations are included in a self-consistent 
manner. In particular, one has to derive self-consistent equations-of-motion for the time-ordered double-time Green's 
function (sub-indexes are omitted here): 



or in terms of their analytic pieces: 



iG(t,t') ee {T c f(t)f(t')), 

iB(t,t') = (T c b(0(t')), (5) 



iG(t,t') = G > (t,t')8(t-t')-G < {t,t')8{t' -t), 

iB(t,t') = B > (t,t')0(t-t') + B < (t,t')0(t' -t); (6) 



A rigorous and well established way to derive these equations-of-motion was first introduced by Kadanoff and 
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Baym [34 j, and has been related to other non-equilibrium methods (like the Keldysh method) by Langreth, see Ref. [35] 
for a review. Here, we just show numerical results of the NCA equations for our problem and refer the interested 
reader to Refs. [3TH331 136] for details. In particular, the density of states is given by 



where A r (e) is the Fourier transform of the retarded Green's function A r (i) = G r ^ La (t)B < (—t) — G< a (t)B a (—t). 
Following Meir and Wingreen in Ref. [37 , the current is given by: 



with A< a (e) the Fourier transform of A< a (t) = iG< a (t)[B r (-t) - B a (-t)} and f a (e) = the Fermi-Dirac 



function at each reservoir held at a chemical potential fi a such that the applied bias voltage is defined as eV = [Ir — {Il- 



In the following, we present results for the density of states and the differential conductance ^ = ^ = —^y- 
Fig. [7] shows the density of states (DOS) for zero bias for symmetric and asymmetric couplings to the reservoirs. 
For B = (zero Zeeman splitting), the peak on the DOS is away from zero frequency while for a finite Zeeman 
splitting B 7^ the central resonance [due to the orbital degree of freedom] shifts to uo — (shown in the Figure for 

Azeeman = 30T^ ^ ') . The results for symmetric and asymmetric couplings are exactly the same because + Tr is 
the same in both cases. 

The behaviour is very different for finite bias. Fig. [8] shows the DOS for B — and different values of the bias 
voltage. First, note that eventually, for a certain value of V, the DOS peak splits in two. V is applied such that the 
left reservoir has hl = V/2 and the right one hr = —V/2 with respect to the Fermi level of the dot. The splitting in 
the DOS is due to the coupling to the two reservoirs displaced with respect to the Fermi energy. The splitting does 
not occur symmetrically because the dot sees the two reservoirs differently, even in the case of symmetric couplings 
Tr = Tl- If a magnetic field is also applied, a Zeeman splitting occurs (see Fig. |9| and, in turn, each new peak splits 
due to the finite bias. 

Let us turn now to the dl/dV curves. For symmetric contacts (T L = Tr = 1) the differential conductance peak 
is always centred at zero bias (see Fig. [lQ^left)). A magnetic field produces a Zeeman splitting and the peak that 
remains at V = is due to the orbital degree of freedom. Remember that the peak of the DOS for the symmetric 
case and B = is away from uj = (Fig. [7] (left)), however, the dl/dV is peaked at zero bias. On the other hand, 
for the asymmetric case (shown in Fig. [To| (right)) the differential conductance peak occurs at finite positive bias for 
B = 0. Applying a magnetic field produces a Zeeman splitting and a shift towards negative values of the bias of the 
central SU(2) peak. The position of this central peak seems to saturate for a sufficiently large value of the magnetic 
field. This can be understood from the density of states: as the magnetic field splits the Kondo resonance, owing to 
the Zeeman splitting of the spin sector, more spectral weight is displaced towards lower energies (close to the Fermi 
energy), see Fig. [9j When the coupling to the leads is asymmetric, this is reflected as a Kondo peak in the differential 
conductance that moves towards lower bias voltages (near zero bias). The peak saturates when spin fluctuations 
are frozen by the magnetic field. This asymmetric coupling situation is the one discussed in the main text to show 
agreement with the experimental data. 
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